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Abstract. We study the structure of colloidal fluids with reference to colloid-polymer 
mixtures. We compare the one component description of the Asakura-Oosawa (AO) 
idealisation of colloid-polymer mixtures with the full two-component model. We also 
\ consider the Morse potential, a variable range interaction, for which the ground state 

clusters are known. Mapping the state points between these systems, we find that the 
pair structure of the full AO model is equally well described by the Morse potential or 
I> ' the one component AO approach. We employ a recently developed method to identify 

in the bulk fluid the ground state clusters relevant to the Morse potential. Surprisingly, 
, when we measure the cluster populations, we find that the Morse fluid is significantly 

closer the full AO fluid than the one component AO description. 



PACS numbers: 82.70.Dd; 82.70.Gg; 64.75. +g; 64.60.My 
1. Introduction 

Although in principle colloidal dispersions are rather complex multicomponent systems, 
the spatial and dynamic asymmetry between the colloidal particles (10 nm-1 /im) and 
smaller molecular and ionic species has led to schemes where the smaller components 
are formally integrated out pQ. This leads to an effective one-component picture, where 
only the effective colloid-colloid interactions need be considered. The behaviour in the 
original complex system may then be faithfully reproduced by appealing to liquid state 
theory [2] and computer simulation [3j. Since the shape of the particles is typically 
spherical, and the effective colloid-colloid interactions may be tuned, it is often possible 
to use models of simple liquids to accurately describe colloidal dispersions. 

Central to this one-component approach is the use of a suitable colloid-colloid 
interaction u(r). Notable early successes include the Derjaguin, Landau, Verwey and 
Overbeek theory of charged colloids [3] and the Asakura-Oosawa (AO) theory of colloids 
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in a solution of polymers 0E], subsequently popularised by Vrij [7j. While theories such 
as these have been used to describe colloidal model systems in which the interactions 
may be tailored with very considerable success [HI El [TO] , the general situation is often 
considerably more complex. 

In the colloid-polymer mixtures of interest here, the effective colloid-colloid 
interactions are set by the polymer chemical potential. One imagines a polymer reservoir 
coupled to the colloidal suspension, in which, if the polymers are ideal as assumed by AO, 
then the polymer chemical potential is proportional to the concentration. In practice, 
experimental systems seldom feature coupled polymer reservoirs, so one is often limited 
to knowledge of the polymer concentration in the sample cell; for a given polymer 
concentration, the chemical potential varies with colloid volume fraction, due to the 
volume excluded to the polymer by the colloids. The volume accessible to polymer is 
also dependent upon phase separation and colloidal crystallisation. In other words, the 
effective colloid-colloid interaction can vary with colloid concentration and also change 
as a function of time, giving rise to novel kinetic pathways and (unlike simple atomic 
substances), a triple coexistence region [11]; meanwhile external fields such as gravity 
may couple with the multi-component nature of the colloid-polymer system to yield 
novel phenomena such as floating colloidal liquids [12]. 

Even in the case of a one-phase colloidal fluid in coexistence with a polymer 
reservoir, for polymer-colloid size ratio q > 0.154 [13], the effective colloid-colloid 
interaction has a many-body component and thus is dependent upon colloid volume 
fraction,while for smaller size ratios the one-component mapping has been shown to 
be exact [Hj. Nevertheless, one may integrate out the polymer degrees of freedom to 
arrive at an effective one-component description for the colloids, as given by AO [6] 
and Vrij [7J. It is worth noting that there is more than one approach to determining 
the effective one-component interaction in a multicomponent system, and that these do 
not always give the same result [15] . The effective one-component description has since 
been extended to include these many-body effects [HJ [T7l [HI [19]. Other important 
departures from the assumptions of Asakura and Oosawa include non-ideal polymer- 
polymer interactions [T5l [20] . which have considerable implications for phase behaviour 
and interfacial properties [21] along with electrostatic interactions between the colloids 
[22]. 

The validity of the one-component approach in describing the colloid-colloid 
interactions has also been investigated experimentally. The interaction between a colloid 
and a glass wall can be accurately measured with total internal reflection microscopy [23J, 
while the interaction between two colloids confined to a line can be measured using 
optical tweezers [231 122] • An alternative approach is to measure correlation functions 
and invert them to extract the effective potentials. Traditionally this has been achieved 
by scattering techniques that measure the reciprocal space structure factor S(k) [21 126]. 
Another means is to determine the structure in real space in 2D and 3D at the single 
particle level using optical microscopy [271 EH], after making some assumptions about 
the system, one may deduce the effective colloid-colloid potential. This may be done 
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with sufficient precision that interaction potentials can be quite accurately determined 
both for purely repulsive systems [28, 29J and for systems with attractive interactions 



The possibility of direct visualisation of colloidal fluids also allows, for example the 
clusters formed to be studied [3U [32} [33] . Lu et. al. explored the idea, introduced by 
Noro and Frenkel in their 'extended law of corresponding states' [34J, that the structure 
of these dilute attractive fluids (the so-called energetic fluid regime [35]) is somewhat 
insensitive to the exact nature of the potential [32]. We have also recently argued 
that the (known) ground state clusters formed by systems interacting under the Morse 
potential (fig. CD) are also relevant to colloid-polymer mixtures. Interestingly, recent 
work suggests that in fact, hard core systems such as colloid-polymer mixtures might 
exhibit somewhat richer (degenerate) topologies of ground state clusters, as more than 
one structure can have identical numbers of bonds [36] . 

Here, we investigate the validity of the one-component approach in colloid-polymer 
mixtures by comparing the full Asakura-Oosawa multi-component model with explicit 
polymers and the one-component AO model [HI [7]. Given a suitable choice of parameters, 
the variable-ranged Morse potential can provide a good approximation to the one- 
component AO potential. In addition to the fact that the ground state clusters are 
known for the Morse potential, we note that its continuous form is amenable to Brownian 
and molecular dynamics computer simulations. We therefore also compare the Morse 
potential by applying the law of corresponding states to map the Morse to the one 
component AO interaction. We consider the structure of the resulting dilute colloidal 
fluids. In addition to conventional pair-correlation function-based methods, we employ a 
recent-developed method which identifies structures topologically equivalent to isolated 
clusters [37J. 

This paper is organised as follows. In section [2] we introduce the simulation 
methodology and our approach for comparing different interaction potentials, our results 
are presented in section [3] and we conclude with a discussion in section [H 

2. Simulations and Interaction Potentials 

The seminal theory of colloid-polymer mixtures is that of Asakura and Oosawa (5[[6]. 
Here colloids are treated as hard spheres with no permitted overlap. Polymers are ideal, 
and may freely overlap with one another, but the polymer-colloid interaction is also 
hard, in that no overlap is permitted. That is to say, the colloid-colloid interaction ucc-, 
colloid-polymer interaction ucp and polymer-polymer interaction upp read 



[30J. 




for r < a 



for r < (cr + o>)/2 
for r > (a + op)/2 



for r > a 



(1) 
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Figure 1. (color online) The ground state clusters for the short ranged Morse potential 
(po = 25.0) for m < 14 particles. Here we follow the nomenclature of Doye et. al. [38] , 

(3u PP {r) = 0. 

where r is the centre to centre separation of the two colloids/polymers and j3 = l/Zc^T, 
where T is temperature, kp is Boltzmann's constant, a and o~p are the diameters of the 
colloids and polymers respectively. 

Some comments on the derivation of the one-component description are in order. 
For a more complete description the reader is referred to Dijkstra et. al. [13]. The 
Hamiltonian of the AO model is thus 

H = Hcc + Hep + Hpp (2) 

where 

No 

H C c = ^2u cc (r) (3) 

N c N P 

H C p = Y J Y, u ^ r ) ( 4 ) 

N P 

H PP = J2upp(r) = (5) 

where Nc and Np are the respective numbers of colloids and polymers. Dijkstra et. al. 
cast the thermodynamic potential F of the colloid-polymer system as 

ex P [-/?F] = Y, TTTTfeTTT / dRN ° / dR N - exp[-^H cc + H p§) 
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N n \A% Nc 



[ dR N ^ex V [-PH EFF ] (7) 
Jv 



where zp is the polymer fugacity Ac, is the thermal De Broglie wavelength of the 
colloids, R Nc and R p are the coordinates of the colloids and polymers respectively. 
H EFF = Nqc + ^ is the effective Hamiltonian of the colloids. 

Now Q is the grand potential of the fluid of ideal polymer coils in an external field 
of N c colloids with coordinates R, and may be expanded as 

tt = tt + fij + Q 2 + ■ ■ ■ (8) 

where Qq is a 0-body term (the Grand potential of an ideal polymer system) Q\ is a 
1-body term related to the volume excluded by the Nc colloids and Q 2 is the two-body 
term. Dijkstra et. al. show that all higher order terms are zero for polymer colloid size 
ratios q = opjo < 0.154 [14] . The two-body term 

n 2 = J2^Ao(r) (9) 

N c 



where 



(3u A o(r) 



-K{2R G fz P {l + qf 



3 



6 q 

for r ><t + {2R G ), 



Now the polymer fugacity zp is equal to the number density ppp of ideal polymers in a 
reservoir at the same chemical potential as the colloid-polymer mixture. Thus within the 
AO model, the effective temperature is inversely proportional to the polymer reservoir 
concentration. The interaction induced by the polymers in equation fflUj) is identical to 
that given by AO [6] and Vrij [Tj. 

We also use the Morse potential which reads 

/3u M (r) = Pe M exp[p (o- - r)]{exp[p (e>" - r)] - 2} (11) 

where po is a range parameter and (3sm is the potential well depth. We set po = 25.0 to 
simulate a system with short-ranged attractions similar to a colloid-polymer mixture. 



2.1. Comparing different systems 

In order to match state points between the Morse and one component Asakura-Oosawa 
interactions, we use the extended law of corresponding states introduced by Noro and 
Frenkel |34j. Specifically, this requires two interactions to have identical well depths and 
reduced second virial coefficients B 2 where 

B* 2 = B 2 /%ia% FF (12) 
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Figure 2. (color online) Interaction potentials used: Morse (blue) and one-component 
Asakura-Oosawa (cyan). Both are scaled by the effective hard sphere diameter (Jeff- 



where cjeff is the effective hard sphere diameter and the second virial coefficient 

oo 

B 2 = 2n J drr 2 [1 - exp (-(3u{r))} . (13) 
o 

The effective hard sphere diameter is defined as 

oo 

oeff = J dr[l- exp {-l3u RE p{r))} (14) 
o 

where the repulsive part of the potential urep is where u(r) > 0. Thus we compare 
different interactions by equating B 2 and o~eff- The latter condition leads to a 
constraint on number density 
_ Nna% FF 

pEFF ~ (15) 

where V is the volume of the simulation box. 



2.2. Simulation Details 

For the one-component systems, we use standard Monte-Carlo (MC) simulation in the 
NVT ensemble j3] with iV = 2048 particles. Each simulation was typically equilibrated 
for 10 7 MC moves and run for a further 10 7 moves. For each state point we performed 
ten independent simulation runs. We confirmed that the system was in equilibrium on 
the simulation timescale by monitoring the potential energy. The Morse potential is 
truncated and shifted at r = 2. ha. In the case of the full AO system, we use Monte- 
Carlo simulation, with polymers included grand-canonically [31 [39]. The interaction 
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potential for the one- component AO is taken as eq. fllUp with the additional hard 
sphere colloid-colloid interaction ucc( r ) [ ec L- ©]• 

We match the Morse and one component AO using eqs. ( fl2l) ( fl5l) by requiring 
the interactions to have the same well depth We set a well depth of 2.0/c#T and colloid 
volume fractions of (pc = iro- 3 pc/6 = 0.05, (pc = 0.25 and (pc = 0.445, where pc is 
the colloid number density. For the Morse interaction, with range parameter p — 25.0, 
this leads to an effective hard sphere diameter o~eff ~ 0.9696a [eq. f fT4|) ]. Applying 
equation ( fl5l) we therefore have a slightly higher volume fraction in the Morse system 
of 4>m ~ 1.10970c- m the one-component AO system, these Morse parameters map via 
eq. ( ITUl) and eq. ( 1121) to a polymer-colloid size ratio of q ~ 0.2575 and polymer reservoir 
number density pp^ « 0.5597<jp 3 . It is worth noting that there is some sensitivity in 
the mapping we have used to the depth of the attractive well. We have taken a value 
of f3eM = 2.0, which we fix throughout this work. However, the 'hardness' of the Morse 
potential depends upon (5em, as, consequently, does the effective hard sphere diameter. 
In principle, one should therefore repeat the mapping for each (3eM- 

The full AO system is challenging to simulate, especially when there is a 
considerable size discrepancy between the colloids and polymers, leading to very large 
numbers of particles in the system [18]. Of course, this is one of the attractions of 
using a one component description. Here we could only equilibrate the system to our 
satisfaction for the higher densities, (pc — 0.25 and (pc = 0.445, owing to the vastly 
reduced number of polymers at higher colloid density. We used N = 256 and N = 512 
for (pc = 0.25 and (pc = 0.445 respectively. The system was equilibrated for 3 x 10 7 
MC moves of either polymer or colloid in each case. Unlike the one component systems, 
two simulations per state point were performed in the case of the full AO system. In 
comparing the full AO system with the one component systems, we only consider the 
colloids and ignore the polymer coordinate data. 

2.3. The topological cluster classification 

To analyse the structure, we identify the bond network using the Voronoi construction. 
Having identified the bond network, we use the Topological Cluster Classification (TCC) 
to determine the nature of the clusters in the bulk fluid [3?]. This analysis identifies 
all the shortest path three, four and five membered rings in the bond network. We use 
the TCC to find clusters which are global energy minima of the Morse potential for 
p = 25.0. These clusters are shown in figure [H We identify all topologically distinct 
Morse Clusters. In addition, for m = 13 clusters we identify the FCC and HCP thirteen 
particle structures in terms of a central particle and its twelve nearest neighbours. We 
illustrate these clusters in fig. [IJ For more details see [37]. We found relatively little 
clustering at the moderate attractions (3e = 2.0 at lower and intermediate densities, 
thus we present TCC results for the highest density studied, (pc = 0.445. 
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Figure 3. (color online) Pair-correlation functions at various densities, (a) Low 
density, <j)c — 0.05, (b) moderate density (<frc = 0.25) and (c) high density {<pc — 
0.445). The g(r) for the Morse potential is shown in dark blue (dashed), the full AO 
system in turquoise and the pale blue is the one-component AO. 



3. Results and Discussion 



We begin our presentation of the results by comparing the pair correlation functions of 
the various systems at differing densities, followed by the TCC analysis. 
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Pair correlation functions are shown in figured At low density, g(r) ~ exp[— /3u(r)]. 
This is illustrated in both cases in figure E(a) (<fic — 0.05), in the form of a strong peak 
at contact, reflecting the short-ranged nature of these attractions. There are some 
minor differences. These are in general consistent with the differences obtained from 
the potentials, figure [2} upon taking the low density limit, g{r) « exp[— (3u{r)]. For 
example, the slighter softer Morse potential leads to a slightly slower decay at r < a. 
Likewise, in the range 1.1a < r < 1.2a, the AO decays to unity rather slower than the 
Morse, reflecting the greater magnitude of the AO in that range. In general, however, 
the agreement between the Morse and AO systems is good. 

We now turn to higher densities, in particular to <pc — 0.25 [fig. [31(b)] . In this 
case, we were able to equilibrate the full AO system in addition to the one-component 
descriptions. Packing leads to a second peak around 2a. Again, we see a similar 
behaviour between the different systems. Significantly, the small differences between 
the g{r)s, comparing Morse to firstly the one component AO and then the full AO, are 
similar. That is to say, the one-component AO, which, for example does not include 
many-body interactions [T6| [T8]. shows discrepancies comparable to the Morse potential 
in its description of the full AO system. 

At the highest density studied (4>c — 0.445), overall we find a similar behaviour, 
as may be seen in figure M^c). This is not altogether surprising, as in dense liquids, the 
structure is well-known to be largely dominated by the hard core [40] . Some differences 
are, however apparent. The Morse system has a weaker first peak, than either the 
one component or full AO systems. This is likely due to the lack of an infinitely hard 
core in the Morse interaction. The first peak notwithstanding, the differences between 
all three systems are comparable. In comparing the one-component AO and full AO, 
our results are compatible with the results of Dijkstra et. ai, who found that g(r)s 
produced from the two descriptions were indistinguishable in the case of q = 0.15 where 
the one-component description is exact [33]. 

We now turn our attention to the cluster populations in the dense system (fig. HJ). In 
all these systems, a range of different clusters are found, with none dominating. Thus we 
argue, that when considering energetically locally favoured structures (i.e. clusters), it 
is important to consider the possibility that more than one topology may be important. 
The overall behaviour between the systems is similar. Among the more populous, smaller 
clusters, the 7A pentagonal bipyramid has a rather low population. However 7A is also 
found as part of larger clusters, notably 8B. According to our counting algorithm, if a 
given particle is part of both a 7A and 8B cluster, it is taken as 8B only. A few particles 
are found as FCC crystal fragments (we found no HCP type environments). 

In comparing these systems we see that the one component AO forms rather fewer 
clusters for 8 < m < 10 than the other systems, and none at higher m. Our statistics 
are necessarily more limited for the full AO system, which we believe restricts our ability 
to determine the population of rarer, higher order clusters. For m < 11, the Morse and 
full AO have rather similar populations, except that the cluster population in slightly 
higher for the Morse system in the case that m > 6. We thus argue that in this respect 
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Figure 4. (color online) (color online) Population of particles in a given cluster, for 
4>c = 0.445. N c is the number of particles in a given cluster, N the total number 
of particles sampled. Here we consider only ground state clusters for the Morse 
po = 25.0 system. Dark blue denotes Morse, turquoise the full AO and light blue 
the one component AO. 

the Morse potential accurately reproduces the full AO model. 
4. Discussion and Conclusions 

We have analysed the pair structure and performed a topological cluster classification 
on a range of model systems for colloid-polymer mixtures. Using the extended law of 
corresponding states [3lj, we have mapped the variable ranged Morse potential to a well- 
known one component model for AO colloid-polymer mixtures. We have also considered 
the full Asakura-Oosawa model. In general, we find good agreement between all three 
systems. The relatively small difference in the pair structure between the slightly soft 
Morse potential and one component AO system seems to be accounted for by noting the 
differences in their functional form (fig. [2]). The small discrepancies exhibited between 
the full AO and the one component systems favour either. That is to say, our g(r) 
results suggest that the Morse potential does as good a job of describing the full AO 
system as the one component AO system. 

Although the pair structure may be very similar between these three systems, the 
topological cluster classification reveals significant differences. In particular, the one 
component AO system forms fewer higher order clusters for m > 8 (8B clusters alone 
account for 20% of the particles in the other systems) and we detect no clusters at all for 
m > 10. In this respect, the Morse potential does a better job than the one component 
Asakura-Oosawa interaction in describing the full AO system. 

Some pointers for further work are considered. Dijkstra et. al. [13, [18] have 
developed an elegant means by which the many-body effects implicit in the full AO 
model are taken into account. It would be most attractive to subject this system to 
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an analysis similar to that presented here. Recalling that we were unable to obtain 
sufficient statistics to calculate a g(r) for the full AO system for <pc — 0.05, we note 
that accelerated MC methods such as the cluster move of Vink and Horbach [39] would 
be most helpful in generating sufficient statistics. 

Moving closer to experiments, non-ideal polymers [20] and electrostatic interactions 
[22] may all impact on these conclusions. We have also considered only a few state 
points. Furthermore, we have neglected polydispersity, omnipresent in experimental 
colloidal systems, which has the potential to alter the results of an analysis similar 
to that carried out here. Coordinate tracking, particularly in 3D experiments based 
around confocal microscopy, is prone to measurement errors of around 0.02 — 0.05<r [30J. 
Work to investigate the sensitivity of this analysis to such experimental considerations 
is in progress. Early indications are that the TCC analysis is surprisingly robust to 
experimental tracking errors and polydispersity. 

The system we have chosen (probably) does not have a stable gas-liquid-coexistence. 
However q = 0.2575 is somewhat above the value of q = 0.154 at which 3-body and 
higher order interactions vanish in the AO model [131 EH]; these effects may be non- 
negligible but the similarity in the correlation functions we measure suggests that the 
effects to not too large, although larger polymers would lead to stronger many-body 
effects [18]. Furthermore, larger polymers lead to such a coexistence between colloidal 
'gas' and 'liquid'. The location of the critical point is known to be strongly dependent 
upon the exact model chosen [39], HI]. Moving closer to the critical point, we expect to 
find different results upon comparing the various models. 

Finally, we have considered equilibrium fluids. The behaviour out of equilibrium 
is most important, particularly in the case of, for example colloidal gels [42]. However, 
we are unaware of suitable simulation models for non-equilibrium studies, except one- 
component descriptions with softened cores [4*3l |4*4"]. and the Morse potential [42J. It 
is almost necessary to use one-component descriptions out of equilibrium, due to the 
degree of computation required. Moreover, Brownian dynamics, appropriate to out- 
of-equilibrium situations, is challenging to implement with hard interactions. Out of 
equilibrium, hydrodynamic interactions may also play a role, and have recently been 
applied to attractive colloidal systems [45] . 
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